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Abstract: The soil respiration rates (Rh) in 6-year-old (young), 17-year-old (middle-age), 31-year-old (mature) Casuarina equisetifolia 
coastal plantations were measured using an LICOR-8100 automated soil CO 2 flux system from May 2006 to April 2007. Results show that 
R\ displayed an obvious seasonal pattern across the observed years. The maximum values of Rh occurred at June and July and the minimum 
at December and January. Soil temperature and soil moisture as well as their interaction had significant effects on the monthly dynamics of 
Rh. The analysis by one-way ANOVA showed that R\ had a significantly exponential relation (p<0.05) to soil temperature at soil depth of 5 
cm, and had a linear relation (p<0.05) to soil water content of the upper 20 cm. The result estimated by the two-factor model shows that soil 
temperature at soil depth of 5 cm and soil moisture at soil depth of 20 cm could explain 68.9%—91.9% of seasonal variations in R h . The or¬ 
der of R\. rates between different stand ages was middle-age plantation>mature plantation>young-age plantation. With the increase of 
growth age of plantation, the Qw of If increased. The contribution of Ri, to total soil surface CO 2 flux was 71.89%, 71.02% and 73.53% for 
the young, middle-age and mature plantation, respectively. It was estimated that the annual CO 2 fluxes from /A were 29.07, 38.964 and 
30.530 t-luf -a" for the young, middle-age and mature plantation, respectively. 
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Introduction 

Soil respiration rates ( R s ), or surface soil CO 2 efflux, is the pri¬ 
mary pathway of soil carbon release to the atmosphere (Raich & 
Tufekcioglu 2000). In recent decades, R s has received more at¬ 
tention because it plays a critical role in global carbon cycle and 
in global warming feedback (Fang & Wang 2007). The R s is 
generally defined as two components, autotrophic respiration ( R . d ) 


Foundation project: The research was supported by Key Projects in the 
National Science & Technology Pillar Program during the Eleventh 
Five-year Plan Period (Nos. 2006BAD03A14-01); Important Science 
& Technology Specific Projects of Fujian province (2006NZ0001-2) 
and supported by the Key Laboratory of south mountain timber culti- 

vation, state forestry administration, P . R. China. _ 

Received: 2009-04-01; Accepted: 2009-05-02 
© Northeast Forestry University and Springer-Verlag 2009 
The online version is available at http://www.springerlink.com 
Biography: XIAO Sheng-sheng (1981 -), Male, Ph.D candidate in In¬ 
stitute of Geographic Sciences and Natural Resources Research, CAS, 
Beijing 100101, P. R. China. (E-mail: xssl9811213@163.com) 
*Corresponding author: Email: yegongfu@126.com 
Responsible editor: Zhu Hong 


and heterotrophic respiration (R h ) (Kuzyakov & Larionova 2005). 
In carbon cycle of terrestrial ecosystem, R h is the main process of 
soil carbon loss to the atmosphere. Thus, quantifying contribu¬ 
tion of i?h to total soil CO 2 effluxes is of great importance and 
prerequisite for estimating the carbon balance in terrestrial eco¬ 
systems (Hanson et al. 2000; Bond-Lamberty et al. 2004). 

Most of the previous studies about R h focus on the inland for¬ 
est ecosystems (Hanson et al. 2000; Ngao et al. 2007; Saurette et 
al. 2008; Gaumont-Guay et al. 2008), and very few concerns 
coastal forest (Tyree et al. 2006). As so far, no available informa¬ 
tion has been found on R h in coastal forest ecosystem in China. 
In coastal regions, forest plantation plays an important role in 
wind combating, sand fixation and local environment improve¬ 
ment. As the main plantation species, Casuarinaceae is widely 
distributed in southern China, especially in Zhejiang, Fujian, and 
Guangdong provinces. It was reported that the area of Casuatitta 
equisetifolia plantations has been up to 3 X 10 5 ha in China 
(Zhong & Zhang 2003). 

In this study, the R h of C. equisetifolia plantation at 3 different 
growth stages was separated and measured in 2006 and 2007. 
The main points of the present study were addressed as: (1) sea¬ 
sonal dynamics of Z? h ; (2) effects of temperature, soil moisture 
and forest age on R h ; (3) contributions of R h to total soil surface 
CO 2 flux; and (4) annual fluxes of i? h in the plantations at three 
different growth stages. 
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Materials and method 

Site description 

The experiment site was situated in Chihu State-owned shelter 
forest of Huian county in Fujian Province, China (latitude 
24 o 30'-24°55'N and longitude 118°40'E—118°55'E,). The study 
site is at elevation of 5-400 m, with an average slope of <10°. 
The climate is maritime monsoon climate, with mean annual air 
temperature of 19.8 °C. The frost-free period is 320 d. Annual 


mean precipitation is 1000-1500 mm, with annual mean evapo¬ 
ration of 2000 mm. The typical soil is homogeneity aeolian sand 
soil with low nutrients, and the layer thickness is 80-100 cm. 
The forest communities are dominated by C. equisetifolici pure 
plantations with few shrubs and grasses. 

Three plantations in different growth stages (6-year-old, P6), 
(17-year-old, PI7), and (31-year-old, P31) were chosen for the 
study. In each plantation, a plot of 20 m*20 m in size was set up. 
Before measurement of soil respiration, the characteristics of 
forest community and soil physiochemical properties were inves¬ 
tigated (Table 1). 


Table 1. The characteristics of forest community and soil physiochemical properties in C. equisetifolia plantations at three different stand ages 




Community 



Soil physiochemical properties (0 

-20 cm) 


Stand ages 

Density 
(trees ha' 1 ) 

Average 

DBH (cm) 

Average tree 

height (m) 

pH 

Soil bulk density 
(g-cm' 3 ) 

SOC (g kg 4 ) 

TN (mg-kg" 1 ) 

C/N 

P6(6-year-old) 

2500 

5.5 

9.5 

5.4740.08 

1.2240.05 

7.1441.14 

72.7243.81 

98.2145.62 

PI 7(17-year-old) 

2202 

12.2 

10.8 

5.2440.05 

1.2440.05 

12.0041.65 

76.9044.12 

156.0547.73 

P31(31-year-old) 

1962 

15.1 

13.7 

5.1040.05 

1.3740.09 

12.4341.69 

81.1544.26 

153.2747.32 


Measurement methods 

In each plantation, the roots were excluded in three subplots (1 
mx 1 m) by using trenching method. Moreover, the total surface 
soil C0 2 effluxes were measured in one subplot (50 cmx50 cm). 

R s and R h were measured at a certain day in the middle of 
every month from May 2006 to April 2007. Measurements were 
at hourly intervals throughout a daytime from 08:00 a.m. to 
18:00 p.m. Every measurement lasted for 2 min with 3 replicates. 
In addition, four observations on daily variations (beginning at 
08:00 in the morning and ending at 08:00 of next morning) were 
carried out for each site on July 2006, October 2006, January 
2007 and April 2007, respectively. 

R s and R h were measured using an LI-8100 Automated Soil 
C0 2 Flux System (LiCor Inc., Lincoln, NE, USA) equipped with 
a short-term chamber. Integrated pump circulates the headspace 
air from the chamber to the non-dispersive infrared (NDIR) gas 
analyzer during the closed state of the chamber and the C0 2 
concentration data as well as the calculated flux rates are re¬ 
corded in the system. A PVC soil collar (inner diameter 20.4 
cm*7.5 cm tall) was inserted into soil at 5-cm depth on May 
2006 in every R s and R h point, and left in the same locations 
throughout the study period. Meanwhile, all living vegetation 
aboveground in the plots was removed (don't disturb the litter 
layer) and the plots were kept free of living plants during the 
study period. 

Simultaneously, air temperature (T a ), surface soil temperature 
(T 0 ), soil temperature at a depth of 5 cm (T s ) and soil water con¬ 
tent at a depth of 20 cm were measured by a DHM2 ventilated 
thermometer, convoluted tubule geo-temperature table, the tem¬ 
perature probe belong to the LI-8100 and drying method, respec¬ 
tively. 

Data analysis 

Daily average flux of R h according to the daily variations in 24 h 

<£) Springer 


was observed in the four days on July 2006, October 2006, Janu¬ 
ary 2007 and April 2007, respectively. It was found that the ef¬ 
fluxes measured during 10:00-12:00 in the morning can basically 
represent for the daily average flux of R h . Kessavalou et al. 
(1998), Du et al. (2006) and Jin et al. (2007) also used similar 
treatments in studying the seasonal dynamics of soil greenhouse 
gases. Therefore, in this study, we used the R h rates during 
10:00-12:00 in the morning to represent for the R h daily average 
flux in the other eight months. 

The relationship between R h and temperature (7) was simu¬ 
lated with an exponential equation and (J 10 value (defined as the 
increase in respiration rate per 10°C increase in temperature). 

R, = ae bT ,Q10 = e' ob (1) 

where, a is the basic respiratory (the respiration rate when tem¬ 
perature at 0 °C), b is the temperature sensitive coefficient. The 
relationship between R h and soil surface water content (W, vol¬ 
ume to volume, 0-20 cm) is simulated with a linear equation. 

R h = cW + d (2) 

where, c and d are the regression coefficients. 

Correlation analysis showed that R h had a better correlation 
with T 5 than that of T a and T 0 . So the dependence of Tj and soil 
water content ( W) on R h was analyzed by a nonlinear regression 
model as follows. 

R h = axe bT xW c (3) 

where, T is the soil temperature at 5-cm soil depth, W is the soil 
water content, and b and c are regression coefficients. 

H c = R„/R s (4) 
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where, H c is the contribution of R j, to total soil surface C0 2 flux. 

Calculating of annual CO 2 flux of R/,: firstly, the average daily 
values of T d (from an automatic weather station of the forest 
center) were transfonned to Tj according to their significant cor¬ 
relations and the average monthly values of Tj were calculated; 
secondly, combinating with the average monthly value W (3 days 
per month for each forest type), Tj was inputted into the nonlin¬ 
ear model (Equation (3)) to get the C0 2 emission of R b in every 
month during the study period, and finally, the fluxes in the 12 
months were cumulated to get the annual flux of Ry,. 

We used SPSS 13.0 (SPSS Inc. 2001) to do the statistical 
analyses. Relationships between R h and T a , T 0 , Tj and W, and the 
differences among 3 plots were analyzed using the One-way 
ANOVA. Two-factor regression analyses were made under the 
nonlinear module. Graphs were prepared by using Excel (Mi¬ 
crosoft Corp. 2003). 

Results 

Seasonal dynamics of R h 

The R h rate in C. equisetifolia plantation had an obvious seasonal 
pattern, with a similar single peak curve for each of three growth 
stages (Fig. 1). In March and April, soil C0 2 efflux stayed at a 
low level, and increased slowly with the increase of microbial 
activity resulted by the risen temperature and the improvement of 
soil water content. It reached the maximum at June (for Pi7) or 
July (for P6 and P31) and decreased gradually with the minimum 
at December (for P6 and PI7) or January (for P31). The coeffi¬ 
cients of seasonal variation were 23.68%, 36.96% and 30.46% 
for P6, P17 and P31, respectively. 

The mean R h rate (pmol C0 2 m V 1 , mean±SE) in the three 


plots was ranked as: P17 (2.76±0.085) > P31 (2.17±0.055) > P6 
(2.00±0.039). Mean multiple comparisons in SPSS, at the 0.05 
significant level, showed that there was a significant difference 
of R h rates between P6 and PI7 (p=0.018), but no significant 
difference was found between PI7 and P31 (p=0.060), and be¬ 
tween P6 and P31 (p=0.593). 



Month 


Fig. 1 Seasonal variations n I' R l: rates for the C. equisetifolia planta¬ 
tions at three different stand ages 

Models of R h to temperature and soil surface water 

In this paper, 7a> T a , T 5 and W were measured and their effects 
on R h were analyzed (Table 2). The influences of temperature on 
R h were different among T d , T 0 and Tj. The R h rates had no sig¬ 
nificant correlations with T. d and T 0 for the three plots except for 
the relationship between T d and R h in the young-age plantation, 
but had significant relationships with Tj in all the three growth 
stages (P<0.05). Tj could account for 51.73%- 62.87% of R h 
variations. 


Table 2. Models of R t] to temperature and surface soil water content for the C. equisetifolia plantations at three different stand ages 




R= ae br 


R=aW+b 


R= 

=a xe^xtR 


Stand age 

Ta 

To 

Tj 




R 2 

R 2 

R 2 

R 2 

a 

b 

c 

R 2 

P6(6-year-old) 

0.5309 

(/>=0.082) 

0.5908 

(/?=0.116) 

0.6287 

(p= 0.000) 

0.5926 (p=0.000) 

1.973 

0.014 

0.114 

0.689 (/>=0.000) 

PI 7(17-year-old) 

0.3239 

(/>=0.002) 

0.4510 

(P=0.067> 

0.5173 

(p=0.000) 

0.7865 (p=0.003) 

3.731 

0.020 

0.422 

0.919(^=0.000) 

P31(31-year-old) 

0.4750 

(p=0.157) 

0.5806 

(p=0.719) 

0.6230 

(p=0.000) 

0.8207 (p=0.001) 

2.262 

0.038 

0.552 

0.901 (p=0.000) 


Notes: all models were tested at the significance level of 0.05 (two-tailed). T a , To. T- indicate air temperature, surface soil temperature, soil temperature at 5-cm 
depth, respectively. W means soil water content at 20-cm depth. The temperature indicator in the two-factor model was based on Tj. 


In addition, the Q\ 0 values of R h based T 5 were 1.39, 1.44, and 
1.95 for P6, P17 and P31, respectively, which showed that the 
temperature sensitivity of microbial respiration improved with 
the increase of stand ages. Kirschbaum (2006) and Fan et al. 
(2008) summarized the influencing factors of R h , and pointed out 
that temperature, moister and the soil substrates were the main 
factors. In our study, no obvious differences between temperature 
and soil moister conditions were found among C. equisetifolia 
plantation at three stand ages; furthermore, soil organic carbon 


contents were similar. Soil heterotrophic respiration is the de¬ 
composition process of soil organic matter and litter dominated 
by soil microbial. Hence we hypothesized that the differences of 
Q i0 values in our study were mostly likely caused by different 
temperature sensitivity of litter decomposition in the plantations 
at different stand ages. Of course, the hypothesis must be tested 
and confirmed in subsequent study. 

The linear models showed that the R h rate had significant cor¬ 
relation with surface soil water content. Soil water content ac- 
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counted for 59.26%-82.07% of R h variations, which stated 
clearly that the effect extent of soil moisture on microbial respi¬ 
ration was larger than that of temperature because of the sandy 
soil of coastal protective forest. 

The biological processes in R s can be considerably influenced 
by soil temperature and moisture, and their interactions corre¬ 
lated to R s must be considered, so do R h (Zhang et al. 2005; 
Wang & Yang 2007), which also was confirmed in our study. 
Soil temperature at soil depth of 5 cm and moisture at soil depth 
of 20 cm had better prediction on R h than using respective inde¬ 
pendent variable, and R 2 values were 0.689, 0.919 and 0.901 for 
P6, PI7 and P31, respectively. 

Contribution of i? h to total surface soil CO 2 efflux and its sea¬ 
sonal pattern 

In P6, the H c varied from 65.50% to 78.56%, with an average of 
71.89%; while in P31, the H c varied from 66.96% to 81.26%, 
with the average of 73.53%. In P17, the H c varied from 64.82% 
to 77.74%, with an average of 71.02% (Fig. 2). However P31 
had higher H c than P6 and PI7 in most months, and there was no 
significant difference of H c among the three plantations at dif¬ 
ferent growth ages at 0.05 significant level (P= 0.326). 

The seasonality of H c showed a similar pattern for the 3 plots 
(Fig. 2). From May 2006 to August 2006, the H c maintained at a 
low level with an average of 67.51%. H c increased and peaked 
on February-March 2007, which showed a reverse seasonality 
pattern comparing to the air temperature. 

In P6, the H c varied from 65.50% to 78.56%, with an average 
of 71.89%; while in P31, the H c varied from 66.96% to 81.26%, 
with the average of 73.53%. In P17, the H c varied from 64.82% 
to 77.74%, with an average of 71.02% (Fig. 2). However P31 
had higher H c than P6 and PI7 in most months, and there was no 
significant difference of H c among the three plantations at dif¬ 
ferent growth ages at 0.05 significant level (P= 0.326). 



Month 


Fig. 2 Seasonal dynamics of the contribution of Rh to total surface 
soil C0 2 efflux in the C. equisetifolia plantations at different growth 

ages (P6 means 6-year-old; P17, 17-year-old; P31, 31-year-old) 

The seasonality of H c showed a similar pattern for the planta¬ 
tions at three different stand ages (Fig. 2). From May 2006 to 
August 2006, the H c maintained at a low level with an average 
of 67.51%. H c increased and peaked on February-March 2007, 

Springer 


which showed a reverse seasonality pattern comparing to the air 
temperature. 

Annual soil CCA flux of i? h 

The results showed that T. d was significantly correlated to T s , and 
the average of R 2 values in the plantations at three different 
growth ages reached 0.92 (p<0.05), so the air temperature could 
be converted to the soil temperature at a soil depth of 5 cm. Ac¬ 
cording to the method mentioned above, we estimated the annual 
soil C0 2 effluxes of R h and the values were 29.072, 38.964 and 
30.530 t-ha' 1 -a‘ 1 for the young, middle age and mature plantations, 
respectively (Table 3). 


Table 3. Linear relationships between J? and T, and annual C0 2 flux 
of R h of the C. equisetifolia plantations at three different stand ages 

(tha'a 4 ) 


Stand age 

Linear relationship between 

T 5 and T a 

Annual C0 2 

flux of Rh 

P6(6-year-old) 

y= 0.927* + 1.724 
(P 2 =0.9333, P<0.05, n=36) 

29.072 

PI 7(17-year-old) 

y = 0.927* + 1.724 
(P 2 =0.9333, P<0.05, n=36) 

38.964 

P31(31 -year-old) 

y = 0.895.v + 2.851 
(R 2 =0.8850, PO.05, n=36) 

30.530 


Discussion 

R h an d the environmental variables 

Our study shows that the mean R h rate in C. equisetifolia planta¬ 
tions was 2.31pmol C02-m' 2 s' 1 , which was higher than that of 
native Castanopsis kawakamii forest (1.88pmol CCA-m^-s" 1 ), 
monoculture Castanopsis kawakamii plantations (0.82|j,mol CO 2 
•m'^s' 1 ) and Chinese fir plantations (0.73|imol C02-m' 2 s' 1 ) in 
Fujian Province (Yang et al. 2006). Zhou et al. (2000) reported 
that the mean R h rate was 0.94(j,mol CCVm^s' 1 in deciduous 
broad-leaved forest and 1.30|amolCC>2'm 2 -s‘’ in evergreen 
broad-leaved forest. 

Many reports about soil respiration pointed out that tempera¬ 
ture and moisture were the main corresponding enviromnental 
factors (Fang & Moncrieff 2001; Rodeghiero & Cescatti 2005). 
R h is also influenced by both climatic factors and human activity 
(Jia et al. 2005). Temperature is the main controlling factor of i? h 
in terrestrial ecosystems and would be regulated by soil moisture 
(Davidson et al. 2000). When soil moisture was the limiting fac¬ 
tor, R h was not controlled by temperature, and instead, soil mois¬ 
ture became the major driver of R h. This phenomenon was con¬ 
firmed by Kucera & Kirkham (1971), Zhou et al. (2004) and 
Raich & Mora (2005). In this study, air and soil temperatures 
were maintained at a relatively high level and had little variation, 
but soil moisture in the coastal sandy soils was very low during 
the study period (averaged 4.33%). Therefore, surface soil mois¬ 
ture was regarded as the limiting factor of R h and presented a 
high relationship with soil heterotrophic respiration. The maxi- 
mums of R h appeared on June or July when there were better soil 
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moisture conditions, compared with soil moisture on August with 
the highest air temperature. 

Contribution of R h to total soil CO 2 flux 

The contributions of R h to total surface soil CO 2 flux varied 
greatly with different plantations and sampling time. Raich & 
Tufekcioglu (2000) revealed that H c had a low value in cold 
biomes, ranging from 7% to 50% in arctic tundra and 11% to 
38% in boreal forests. In temperate regions, the H c ranged from 
50% to 67% in broad-leaved forests from 38% to 65% in pine 
forests. Therefore, the trend of increase of H c with temperature 
increasing in macro-scale is quite clearly. The annual mean air 
temperature is 19.8°C at this field, so the average H c at 72.15% 
in the three C. equisetifolia plantations was reasonable. 

However trenching method was a mature approach to separate 
R a and R h, which also could affect R h and its contribution to total 
surface soil CO 2 flux; for example, Uchida et al. (1998) had 
pointed that, in the trenched plots, the decomposition of the root 
residuals could largely enhance the soil heterotrophic respiration 
and elevate the contribution of soil microbial respiration to total 
soil CO 2 effluxes. Moreover, temperature and soil moisture 
would change in the trenched plots, which could also alter the 
soil microbial activity (Wang & Yang 2007). So far, however, the 
estimated range of trench effect and how to reduce or even 
eliminate the trench effect have been not well studied. 

In the C. equisetifolia plantations, the H c varied at different 
growth seasons. In winter and spring, H c presented a high value; 
while in summer and autumn, H c was decreased. The results 
were agreed with the findings by Zhang XQ et al. (2005) and 
Chang et al. (2007). Li et al. (2000) reported that soil respiration 
was mainly contributed by the decomposition of soil organic 
matter when air temperature was below 15°C; but when air tem¬ 
perature was above 15°C, the contribution of plant root respira¬ 
tion was larger than microbial respiration. In our study, the H c 
decreased when air temperature and soil temperature increased. 
This phenomenon may be attributed to the booming of root 
metabolic activity and the higher (J 10 of the root respiration than 
that of microbial respiration (Boone et al. 1998; Epron et al. 
2001). In addition, the higher H c in winter and spring was partly 
due to the frequent windy weathers which increased the litter of 
branches. 

Effects of stands at different ages on R h 

Previous studies about the influence factors of R h in forest com¬ 
munities focused upon soil temperature and moisture (Wang & 
Yang, 2007; Bauer et al. 2008; Gaumont-Guay et al. 2008; Graf 
et al. 2008), the extent and mechanism of the effect of stand age 
on R s and R h were not well known. Saurette et al. (2008) reported 
that there was no significant difference of R h between plantations 
at different growth ages. However, the significant difference of 
R h was found in this study (p<0.05), the R h was highest in the 
middle-age plantation and then followed by the mature plantation. 
The lowest R h rate occurred at the young-age plantation. Differ¬ 
ent growth ages of plantation can influence R s and Rh by the type 


and number of microbial communities (Gao et al. 2007), above¬ 
ground biomass (Chapin & Ruess 2001), Quantity and quality of 
litters (Liu & Fang 1997) and soil quality (Shao et al. 2005). The 
differences of Rh rates among different forest communities that 
had similar climatic conditions and soil environment were 
mainly driven by different quantities and quality of soil respira¬ 
tory substrates (Ma et al. 2003; Yang et al. 2005). In our study, 
the C. equisetifolia plantations at three stand ages have the same 
meteorological conditions and soil physiochemical properties; 
especially there were no obvious differences of soil organic car¬ 
bon contents which had low background levels and low contribu¬ 
tion of soil mineral respiration to soil heterotrophic respiration. 
Moreover, the litter quality among the three plantations was 
similar. Thus, we hypothesize that the differences of Rh rates 
between the plantations at three growth ages may be attributed to 
the different litterfall quantity (especially different quantities of 
fresh litter). This hypothesis need proof in the future study. 
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